Exactly solved models of interacting dark matter and dark energy 
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We introduce an effective one-fluid description of the interacting dark sector in a spatially flat 
Friedmann-Robertson- Walker space-time and investigate the stability of the power-law solutions. 
We find the "source equation" for the total energy density and determine the energy density of 
each dark component. We study linear and nonlinear interactions which depend on the dark matter 
and dark energy densities, their first derivatives, the total energy density with its derivatives up 
to second order and the scale factor. We solve the evolution equations of the dark components for 
both interactions, examine exhaustively several examples and show cases where the problem of the 
coincidence is alleviated. We show that a generic nonlinear interaction gives rise to the "relaxed 
Chaplygin gas model" whose effective equation of state includes the variable modified Chaplygin 
gas model while some others nonlinear interactions yield de Sitter and power-law scenarios. 
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I. INTRODUCTION 

At large scales, there are strong evidences for a spa- 
tially flat and accelerating universe transiting from a sce- 
nario dominated by matter accumulated by purely at- 
tractive gravitational effects to another dominated by a 
dark energy component dispersed by repellent gravita- 
tional effects and characterized by a negative pressure 
This behavior has lasted till present and most likely 
continue for ever. The reasons, causes and details of when 
this transition happened, it was not still understood. For 
a sufficiently intense acceleration one can speak rather 
of superacceleration, then the possibility exists that the 
universe has a catastrophic end with a suddenly future 
singularity at a finite time (Big Rip) and a total disinte- 
gration of the well-known structures. 

Astronomical observations suggest that the Universe 
entered an accelerated expansion stage when the value of 
its scale factor was approximately one-half of the current 
one. This important discovery, which was based on the 
observations of the brightness of a class of supernovae 
(SNIa) 0, has been confirmed by precise measurements 
of the spectrum of the cosmic microwave background 
(CMB) anisotropics [|[ as well as the baryon acoustic os- 
cillations in the Sloan digital sky survey luminous galaxy 
sample Q ■ The consensus between cosmologist points in 
the direction that the understanding of the phenomenon 
will probably require a unified comprehension of the grav- 
itational and the other fundamental interactions. 

We will consider fundamentally interacting and uni- 
fied models to investigate the mechanisms governing the 
dynamics of the evolution of the universe from its early 
stage until its recent accelerated phase. 

In the interacting models, the source of Einstein equa- 



tions includes an aggregate of different material fluids 
that interact among them. This is in principle the sim- 
plest, and perhaps the most obvious hypothesis, and it 
is in fact, the one that has provided more advances in 
the knowledge of the phenomenon of the recent accel- 
eration of the universe @, and references therein. 
Observational evidences show three fundamental compo- 
nents: baryons, dark matter and dark energy. Given the 
dynamical similarity between baryons and dark matter, 
it is possible to make a simplified model replacing both 
components with a nearly prcssurclcss dust. 

In unified models, the Einstein equations will have a 
single component working as dark matter and dark en- 
ergy at different stages. It interpolates smoothly between 
a matter dominated phase in the early stage and dark 
energy in the late stage of the evolution, so inducing an 
accelerated expansion of the universe. The Chaplygin 
gas and its extensions were the unified models that have 
been more studied in the literature 0]~ E3 - More re- 
cently, several generalizations of those models have been 
introduced, as for instance, the variable modified Chap- 
lygin gas model with an equation of state, depending 
explicitly on the scale factor [lj|]- 23 1. 
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In the following we will use either interacting or unified 
models indistinctly because there is some sort of resem- 
blance between these different models. 

In section II we consider a dark sector model with en- 
ergy transfer, develop an effective one-fluid description, 
find the source equation for the total energy density and 
investigate the stability of the power-law solutions. In 
section III we introduce the "linear interaction" , impose 
the stability condition, describe some simple interacting 
models and find the exact scale factor as well as the effec- 
tive barotropic index. Also, we present a "general linear 
interaction" which induces a generalized ACDM model. 
In section IV we investigate a "nonlinear interaction" , 
solve the source equation and show that the effective 
equation of state of the dark sector includes several gen- 
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cralizations of the Chaplygin gas. Finally, in section V 
the conclusions are stated. 



II. DARK SECTOR EVOLUTION 



dynamics of p c and p x is dictated by the scale factor in- 
stead of H. Differentiating the first or the second Eq. ([5]) 
and combining with the Eq. ^ or with the Eq. ([7]), we 
obtain a second order differential equation for the total 
energy density 



A. Effective one-fluid description 

We consider an expanding universe modeled by a mix- 
ture of interacting dark matter and dark energy with 
energy densities p c and p x , and pressures p c and p x re- 
spectively. Due to the energy transfer between both dark 
components, they do not evolve separately Q and the 
Einstein equations read 



3H~ 



Pc + Px, 



Pc + Px + 3H(p c +p c 



■P*) = 0. 



(1) 



(2) 



where a is the scale factor and H = a/ a. The conserva- 
tion equation ([2j evidences the interaction among the 
components admitting the mutual exchange of energy 
and momentum. 

For the two dark components we assume equations of 
state p c = (7 C - l)p c and p x = (7^ - l)p x , where the 
barotropic indices j c and "f x are constants. The dark 
matter is composed of nearly prcssurcless components 
with a barotropic index 7 C «1 and the dark energy vio- 
lates the strong energy condition, p x + 3p x < 0, so that 
7a; < 2/3 < 7 C . The total energy density p and the con- 
servation equation for the interacting two-fluid model are 



P = Pc + Px, 



-IcPc -IxPx, 



(3) 



(4) 



where the prime indicates differentiation with respect 
to the new time variable ' = d/dr) = d/3Hdt = 
d/d In (a/ ao) 3 and ao is some value of reference for the 
scale factor. Solving the system of equations (O-O we 
get p c and p x as functions of p and its derivative p' 



Pc 



IxP- 



A 



Px 



IcP + P 



(5) 



where A = 7 C — "f x > is the determinant of the linear 
equation system ©-(U]). 

At this point, we introduce an energy transfer between 
the two fluids by separating the conservation equation 
dH) into the two equations 



Pc + IcPc 



P x + IxPx = Q- 



(6) 



(7) 



Here, we have consider a coupling with a factorized H 
dependence 3HQ where the interaction term Q gener- 
ates the energy transfer between the two fluids. So the 



p" + (ic + lx)p' + JcJxP = QA, 



(8) 



that we call the "source equation" [241 ]. 

The interacting two-fluid model has been reduced to 
an effective one-fluid model with total energy density p 
and total pressure p = p c + p x , whose effective equation 
of state is 



Pip, p') 



-P-P 



(9) 



This enable us to assume an effective one-fluid descrip- 
tion with equation of state p — (7 — l)p and effective 
conservation equation p' + ~/p = 0, where the effective 
barotropic index 7 = (j c Pc + IxPx)/ P ranges between 

lx < 7 < 7c- 

Resuming, given an interaction Q, the total energy 
density p of the effective one-fluid model is determined 
by solving the source equation Eq. (|8|). Once we know 
p, we are able to find the effective equation of state from 
Eq. © and the scale factor by integrating the Friedmann 
equation 3H 2 = p, without knowing p c and p x separately. 
Both energy densities are easily calculated by replacing 
p and p' into the Eq. ([5]). For instance, in the no in- 
teraction case, Q — 0, the energy density of the effective 
one-fluid model is p — c /a 37c + b / a 37x . For any value 
of the constants c and b, p —> b /a 37x , the scale factor 
a — > f 2 / 37 * and the power-law solution i 2 / 37 * becomes an 
attractor. Throughout the paper Cj, bi, ji, ki and rii with 
i = 1,2,3...., will represent constants while c and b will 
be integration constants. 

Basically, we have shown that an interacting two-fluid 
model can be seen as an effective one-fluid model or 
equivalent ly considered as a unified one. Its dynamics 
is given by the two independent Einstein equations 



3H 2 = p, 



3H(p + p) = 0. 



(10) 



These equations cannot determine the three quantities a, 
p, and p because there is one degree of freedom. Usually, 
the system of equations (|10p is closed with an equation 
of state p = p(p). On the other hand, when we split the 
effective energy-momentum tensor into two dark compo- 
nents, the Eqs. (|10j) become Eqs. |T]) and ©, which 
cannot determine the five quantities a, p c , p x , p c and p x . 
Thus, we need to introduce an equation of state for each 
dark component p c = (7 C — l)p c and p x — (j x — l)p x to 
preserve the one degree of freedom of the unified model. 
Then, by changing these equations of state we obtain a 
very large set of interacting models which are equivalent 
to a unified one, meaning that the decomposition into 
dark matter and dark energy is not unique. 
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B. Asymptotic stability 



Hence, when the function F fulfills the two conditions 



The knowledge of stable power- law solutions a = t 2 / 3js 
for the interacting two-fluid model is very useful because 
these attractor solutions determine the asymptotic be- 
havior of the effective barotropic index 7 = —2H/3H 2 = 
7 S . These solutions are associated to a universe ap- 
proaching to a stationary stage i.e., existence of attractor 
solutions 7s . Consistently 7 tends to the asymptotic con- 
stant value 7 S . Then on the attractor 



7s 



IcPc 



IxPx 



Tslc + lx 



(11) 



we conclude that the ratio r = p c j p x turns asymptoti- 
cally constant r s — p cs /p xs , alleviating the problem of 
the coincidence: are the proportions of matter accumu- 
lated by gravitativc effects and dark energy comparable 
at the present time for a strange coincidence or for a fun- 
damental reason? [24[-[29||. For the case 7 S =0, we have 
a final de Sitter stage, H = const, with r s = —i x /lc- 

To investigate the stability of the constant solution y s 
we use the evolution equation of 7, which is obtained by 
replacing p' = —jp and p" = (j 2 — r y')p into the source 
equation ([5]). so we have 



i - (7 - 7c) (7 - ix) = -— 

p 



(12) 



We proceed assuming that a constant solution 7 = 7 S 
of the Eq. (|12|) with -f x < j s < j c exists, and after that, 
we will impose the stability condition so that 7 S be sta- 
ble. An interaction satisfying the existence requirement 
belongs to the class 



Q(7s) = A x ( 7 s - 7 c )(7s - lx)p, 



(13) 



with Q(j s ) < 0. The negative value of Q(i s ) on the 
attractor indicates that the energy is being transferred 
from dark energy to dark matter. It assures that power- 
law solutions a = i 2 / 3 ^ , obtained after integrate 7 S = 
—2H/3H 2 , are stable and the ratio r tends asymptoti- 
cally to the constant value r s . 

The analysis of stability will be performed for interac- 
tions that have the form Q = Q(p c , p x , p' c , p' x , p, p' , p"). 
By using the Eq. (O and p' = —jp with p" = (j 2 — y')p, 
we obtain that p c>x = p C)X (p,p') and p' c x = p' CiX (p' , p"), 
then the interaction becomes Q = Q(i,i', p). For 
simplicity we adopt separability of Q, that is, Q = 
Q(l,l',p) = pQ(l,l') and write 

Q( 7 , 7', P) = A- : ( 7 - 7c)(7 - lx) F( 7) i')p, (14) 

where the function F depends on 7 and 7'. Several in- 
teracting models analyzed in the literature are described 
by the interaction term (|14p . see for instance (2rl|-[37|. 
Other kinds of interactions were explored in [38| . 

By combining Eqs. (|T2|) and (|14[) we rewrite the evo- 
lution equation of the effective barotropic index as 



and 



0) 



F(7 = 7s ,7 / = 0) = 1, 



(is - 7c) (is - 7 a )F 7 (7s,0) 
1 + (is - lc)(ls - lx)F Y (i s ,0) 



(16) 



< 0, 



(17) 

where F 7 and Fy stand for the partial derivatives of F 
with respect to 7 and 7' respectively, then 7 S is a stable 
solution or an attractor. In other words, when the con- 
dition (fTrl)) is satisfied y s becomes a constant stationary 
solution of the Eq. (jT5j) and it is stable whenever the 
stability condition (|17[) is fulfilled. 



III. LINEAR INTERACTION 

Interaction terms depending linearly on the energy 
densities p c , p x and the total energy density p have been 
investigated in a broad class of cosmological models [2~ij - 
[35| . When the dark matter and dark energy are coupled 
with p x , p or linear combinations of p c and p x , the ratio 
r may tend asymptotically to a constant value. This at- 
tractive feature encourage us to explore an extended set 
of interactions. From Eq. ([5]), we have seen that p c and 
p x are linear functions of p and its derivative p' . Tak- 
ing into account that the latter term includes the first 
derivative of the energy densities p' c , p' x , and from Eq. 
([5]), these terms introduce a dependence with p" , we will 
proceed to study the case of considering the "linear in- 
teraction" (Ql) which admits new terms proportionals 
to p' c , p' x and p" . Then, we start building the "linear 
interaction", Ql, as the following linear combination 

Ql = cip c + c 2 p x + c 3 p' c + c 4: p x + c 5 p + c 6 p' + c 7 p", (18) 



in such a way that it verifies the conditions (p~6]l - (|17[) . 
The interaction Ql is also motivated for the fact that 
the source equation ([8]) becomes a linear second order 
differential equation for p. As far as we know, this this 
kind of coupling was not investigated in the literature, 
so we will analyze in detail the Q l . The particular cases 
Qt, = ci Pc + C2p x were previously 

3C|, jag with Q L = c 5 p in Mj-p 

34 1 and with Ql = c-ipx in [26|-[29|. 

By using the Eqs. (HJ) in the Eq. (fl"8)) . the Ql can be 
reduced to a linear combination of the basis elements p, 
p' and p". Finally we get 



investigated in Refs. 
I, with Ql = cip c in 



hp + hp' + hp", 



(19) 



where the constants h are linear combinations of the con- 
stants c, in the Ql (fT5|) . 



h = A 1 (-i x ci + 7 c c 2 + c 5 ), 



(20) 



i = -(l-lc)(l-lx)(F-l). 



(15) 



A 1 (-ci + c 2 - i x c 3 + 7 c c 4 + c 6 ), 



(21) 
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b 3 = A- 1 (-c 3 + C4 + c 7 ). (22) 

Combining p' = — 7 p and p" = ( 7 2 — l')p with the Eq. 
(|T9| . and comparing the Q L ((HI) with the Eq. (JTU we 
obtain the following function F( 7 , 7 ') 



F 



{l - 7c) (7 - 7a) 



[61 -6 27 + 6 3 (7 2 -7')] • (23) 



The existence of a stationary solution 7s of the Eq. (|15[) 
is linked to that the function F satisfies the condition 
(UHJ) leading to the constrain 



A 



{is -7c) (is ~lx) 



\bi - b-2ls + &37.s 



1. 



(24) 



for the coefficient b\, 62 and 63. By solving Eqs. (|24[) for 
62 and inserting it into Eq. (jTTJ)) , we obtain the final form 
of the Q L 



- (7s lc)[ J S lx) p'+hp". 

(25) 

By inserting this Ql into the Eq. (fT2|) we find the two 
constant solutions 



A 



7 



7s 



7 



7e7a - biA 
ls(l-b 3 A) : 

while the stability condition ([TT|) gives 

7 S - 7+ < 0, 



(26) 



(27) 



so, the constant solution 7s is asymptotically stable pro- 
vided that 7s < 7 + with the additional requirement 
lx < Is < 7 + < 7c, yielding a stable cosmological model 
with the asymptotic power-law expansion a = i 2 / 37s . 

The exact general solution of the source equation (|8]) 
for Q l (f2"5"j) and the effective equation of state © are 



p L = c«- 37 -+i»«" 37+ , 



(28) 



universe begins with a mix of dark matter (|30[) and dark 
energy (|3f [) represented approximately by the unstable 
p c oc ( 7 + — 7; C )a~ 3 ' 7+ and p 2 - cx ( 7c — 7 + )a~ 37+ respec- 
tively. After that, the instability of the constant solu- 
tion 7 + induces the universe to evolve from that unsta- 
ble era, characterized by r + = ( 7 + — r y x )/(j c — 7 + ), to 
a stable final stage where the dark matter and dark en- 
ergy densities are dominated by the stable components 
p c oc ( 7s — 7:z)a -37s and p x oc ( 7c — 7s )a~ 37s . The stable 
solution "f s is associated to an asymptotically stable ra- 
tio r s = ( 7s - lx)/{lc - Is) with r + > r Sl showing that 
the linear interaction alleviates the problem of the coin- 
cidence. In turn, the scale factor interpolates between 
the unstable stage, evolving as a oc t 2 / 37 , and the sta- 
ble stage evolving as a oc £ 2 / 37s . Meanwhile the effective 
equation of state (|2T)|) plays the role of a peculiar fluid. 



A. Linear examples 

We will study simple examples where the dark compo- 
nents interact with each other successively by considering 
separately only the terms p c , p x , p and p' of the Ql ([18]) 
and review some of the models investigated with these 
particular couplings. Also, we select the four function 
F(j) in such a way that the condition (fTB)) is satisfied 
identically. 

(1) Cl ^ 



{is ~ 7c)Pc 



F„ 



Is - 7c 

7 -7c 



(32) 



with Q Pa < [341 ]. The stability condition ([i"7]) . is not 
satisfied and the power-law solution a = t 2 / 37s becomes 
unstable. This model contains serious instabilities on the 
perturbations of the dark energy component [35j . 
In this case the solution of the source equation © is 



P P c 



caT 3 " 1 ' + ba~ 3 ^ 



(33) 



PL = {is - 1)pl + (7 + - ls)ba 



-3^ 



(29) 



The interacting model is finally realized when the general 
solution (f2"5)) is inserted into the energy density of each 
dark component ([5]) 

{ ls - lx )ca~ 3 ^ + { 1 + - lx )ba~ 3 ^ + ^ 
Pc = 7 , (30) 



(7c-7s)ca- 3 ^ + (7c~7 + )fra- 37 
Px = ^ , (31) 

For large scale factors the quantities p, p Cl p Xl p' c , p' x , 
p' , p" , p and the Ql behave as a _37s while in the initial 
regimen the above quantities behave as a~ 3j . 

The coupling between the two dark components mod- 
ifies typical characteristics of p c and p x . In fact, the 



For any value of the initial conditions c,b and for large 
scale factor, the total energy density (I3"3l has the limit 



Ppc 



,37, 



meaning that a 



2/3 7x 



since j x < 7 S . 



The model seems to be completely dominated by the in- 
put source p x . In fact, the energy densities ([5]) arc 



Pc 



= A- 1 {is -la,) co"**, 



p x = A 1 (7 C - is) c a 



-3 7s 



-37 x 



(34) 



(35) 



and the ratio r„ 



oc a 



-3(7s-7x) 



— > 0. Then, at late times, 



the interacting two-fluid model with energy transfer Q Pa 
is not satisfactory. However, this coupling can work when 
it is combined linearly with some of the other parts of Ql ■ 
(2) c 2 



-{is ~ lx)Px, 



F„ 



ls-lx 
1-lx' 



(36) 
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with Q„. < 0. This interaction was examined in several 
papers j2r|-j29j. Now the stability condition (JTTJ) is sat- 
isfied and the solution -f s is stable. By solving the source 
equation ([5]) for Q Pa , , we obtain 



~ 3 ^ +ba- 3 ^. 



(37) 



For any value of the initial conditions c,6 and large scale 
factor, the total energy density p — > ci/a 37s and a — > 
t 2 / 3 7s because 7 S < 7 C showing that the interacting model 
is dominated by the attractor 7 S . 

The dark matter and dark energy densities ([5]) are 



Pc 



A" 1 (7s 



Px 



A- 1 ( 7c 



7« 



-3 7s 



(38) 



(39) 



showing that the ratio r p 



r s = (7s - lx)/ (7c - 7s) on 



the attractor. Then, the interaction Q Pm may represent 
adequately an interacting dark sector model which can 
be adapted to the observations. A cosmological model 
with the above Q Psc was proposed for the current uni- 
verse which consists of noninteracting baryonic matter 
and interacting dark components [26| . The evolution of 
a viscous cosmology model was also analyzed by employ- 
ing the energy transfer Q Px [28|. 



(3) c 5 ± 



(7s -7c) (7s -7x) 



F n 



(7s ~7c) (7s -lx) 
il - lc)(l ~ lx) 



(40) 

with Q p < 0. It produces a transition from a dark mat- 
ter dominated phase to an accelerated expansion phase 
dominated by dark energy [3lT|- [33| . 

By imposing the stability condition (|17[) on the func- 
tion F p (|40|) we find that 7 = 7 S is an attractor provided 
lx < 7s < (7c + 7z)/2 < 7c- Inserting Q p into the source 
equation ((8j) we obtain the total energy density 



Pp 



co -37 ' +ba- 3 ^" + ^-" /s) . 



(41) 



For any value of the initial conditions c,b and large scale 
factor, p — > c/a 3ls and the power-law expansion a — > 
t 2 / 37s becomes asymptotically stable. 

The dark matter and dark energy densities (0 arc 



(7s - lx) ca-*i- + ( 7 c - 7s) 6a- 3 ^+T«-T-) 



A 



px 



( 7c - 7s) ca~ 3 ^ + ( 7s - lx ) 6a-3(7c+7x-7.) 



(42) 



(43) 



Thus, the ratio r p tends to r s = (i 3 — l x )/(lc~ls)i being 
r s an attractor. 
(4) c 6 ? 

n (7c - 7s)(7s - lx) , „ l(ls ~7c) (7* -7x) 
Qp' = — r P, F p , - 



7s A 



7s(7-7c)(7~7x) ' 
(44) 



From p' = — 7p < 0, we see that Q p > < is negative. 
Besides, by imposing the stability condition (JTTJ) on the 
function F p i , we obtain that 7 S is an attractor provided 
Is < Iclx- Solving the source equation © for Q p i, we 
find the total energy density 



Pp' 



ca 



' 3 ^ + ba~ 3 ^ h \ 



(45) 



Then, whatever be the initial conditions c and b the total 
energy density has the limit p —> c/a 3ls for large scale 
factor, evidencing that 7 S is an attractor. 

The dark matter and dark energy densities ([5]) are 
given by 



(7s - lx) ca~ 3 ^ + 7.77 1 ( 7c - ls )ba~ 3 ^h, 



A 



(46) 



Px = 



(7 C - 7s) ca" 3 ^ + 7C7J 1 (7s ~ 7.) ba-^'/f' 



A 



(47) 

and r p i = (is — lx)/(lc — 7s) on the attractor. As far 
as we know, this interacting model was not investigated 
in the literature. It appears as a feasible candidate to 
be considered for to describe the evolution of the dark 
sector. 

Finally, by solving the Friedmann equation 3H 2 = pl 
for the source (|28|) . we find the exact scalar factor 



cll = [oj sinh At] 



2/3(7 + -7s) 



t = 



V3b{i+-is) 



' sinh At 



7s/ (7 -7s) 



(48) 



dr, (49) 



where uJ 2 = c/b, see Ref. (39|. Due to 7 + > 7 S , the latter 
equation shows that the variables t and r have the same 
asymptotic limits. Then, we use r instead of t to analyze 
the scale factor and the effective barotropic index 



1L 



7 



7 S sinh ojAt 



cosh wAt 



(50) 



in the two asymptotic regimes. So as t grows the model 
interpolates between the initial i + and the final i s val- 
ues. Eqs. (j4*g]) -([5T) |) allows us to express the total energy 
density (|2"5|) , the dark matter and dark energy densities 
([S]), the ratio tl and the effective pressure © as func- 
tions of the new time r. In particular, at early and later 
times, the asymptotic limits of the ratio 77, become 



T 



lx 



lc-T 



Is - lx 

7c - 7s 



(51) 



which satisfy the crucial relation > r s . 

B. General linear interaction and ACDM model 

We complete the subject of linear interaction by en- 
larging the basis elements with a constant, so that, the 



G 



new base will be Qq/A, p, p' , p" , with Qq/A a constant. 
Although, the effective one-fluid model is able to mimic 
the essential features of a de Sitter scenario, clearly, the 
introduction of interactions could produce different alter- 
natives to this scenario. 

We introduce the "general linear interaction", 



QgL — + Ql, 



(52) 



where the constrain (|24p holds for the constants of Ql- 
A particular type of the general linear combination (|52[) , 
Q = Cq + C\p c + C 2 p x , was analyzed in [3^|. A more 
general linear interaction can be built including a f{rf) 
term to the Q 9 l (|52p . sec Rcf. [24j |. A particular case 
where the interaction Q 9 l oc f(rf) was examined in [40j |. 
Combining Eqs. ©, ^ and ([55]) the Eq. © 

becomes 



(1 - hA)p" 



Is 



Iclx 



(b 1 + b 3l 2 )A]p' 



hdx - hA)p = Q Q 



whose general solution is 



PgL 



A 



iff 



-3 7a 



ba 



-37+ 



(53) 



(54) 



where A e ff — Qq/(j c Jx — b\A) > is the effective cos- 
mological constant, induced by the constant term Qq/A 
in the Q 9 l ■ 

At late times the total energy density has the limit 
p g L — > A e y / and the effective equation of state becomes 
p s=! —A e ff. Thus, the effective one-fluid model can be 
associated with a unified dark sector model whose scale 
factor interpolates between a power-law phase and a dc 
Sitter stage H = Hq = y/3/A e ff being Hq an attractor. 

From Eq. (|5|), the dark energy densities are 



PcgL 



IxA 



eff 



Pc, 



PxgL 



7cA 



eff 



(55) 



where p c and p x are given by Eqs. (f3Dj) and (pTTI) . When 
the total energy density tends to A e // for a — > oo, the 
dark matter energy density has the final limit p cg L — s- 
—jxAeff/A < 0. To relieve this problem we may assume 
a phantom equation state for the dark energy, with y x < 
0, so the energy densities (f5"5"j) become positive and the 
ratio r g L tends to r w = ~^ x /lc > 0. 



IV. NONLINEAR INTERACTION 

Let us assume that the energy transfer between the 
dark matter and dark energy components is produced by 
the following "nonlinear interaction" , 



jiPl 



j2PcPx 



hpl 



f(v)p p 



(56) 



where the Ql is giving by Eq. (fT^j) . /(??)p I/ is a nonlin- 
ear atypical term proportional to a well-behaved function 



/(ry), which depends on the scale factor 7] = In a 3 , with v 
a constant. 

By using Eq. ([5]) we see that the three terms in the 
numerator of the nonlinear part of the Eq. (|56[) , become 
a linear combination of p 2 , pp' and p' 2 . Rearranging all 
these terms, we obtain the final form of the Q n L (l56l) . 

n _ k lP ' 2 + p(k 2 p + k 3 p' + k iP ") + /(7 ? K +1 

where the fcj are combinations of the constants in Eq. 
([56"]) . From Eqs. ([5]) and (|57|) . we obtain the nonlinear 
differential equation for p 



PP 



7c + lx - k 3 , 

— i — ; pp 



l-k 



Iclx - k2^ 2 = f 



1 



1 



(58) 



In the particular case where k^ = 1 and / = 0, the first 
and last terms in the RHS of the Q u l (|5T[) vanish si- 
multaneously, Eq. (|58[) becomes a homogeneous linear 
differential equation for p. However, in other cases, the 
general solution of the source equation (|5"5|) will be ob- 
tained from a nonlinear superposition of the two basis 
solutions. 

Renaming the four constants in the Eq. (|58p by n\, 
7T-3 and ri4 respectively (from left to right) and changing 
to the new variable x = p( 1+n2 ) with n 2 ^ — 1 or ki+k± ^ 
1, the Eq. ([5"5j) turns into 

(59) 



v" + ni .x' + n 3 (l + n 2 )x = n 4 (l + n 2 )f(r]), 



where we have chosen v = — n 2 = fci/(l — k±). If X\h 
and x 2 h are the two basis solutions of the homogeneous 
Eq. (f5T)|) . then the general solutions of the Eq. (|5"5|) can 
be written as a nonlinear superposition of these basis 
solutions p n L = {cxih + bx 2 h + x p ) 1 ^ 1+,l2 \ where x v is 
a particular solution of Eq. ([511]) . From now on, the Eq. 
([5T)]) will substitute the source equation §E§ for the Q n L- 



A. The homogeneous case 



In the / = case, the general solution of Eq. (|59l) is 



Xh = ca 

variable ph 
we find 



ba 



3 A" 1 



Ph 



l/(l+n 2 ) 



Then, coming back to the original 
and using the Eqs. ([5]) and (O, 



,3A~ 



b a 



.'iA 1 



+ 11/(1+112) 



Pch = -D 



Pxh = D 



-,3A^ 



[A +7 a (l +n 2 )]p + bAX — 

p™ 2 



[A +7 c (l + n 2 )]p + 6AA 



,3A^ 



(60) 



(61) 



(62) 



7 



Ph 



1 



A" 



Ph 



bAX a 3 



(63) 



where A - , A + are the characteristic roots of the Eq. ([55)1 



A T = 



-ni T \Jn\ - 4n 3 (l + n 2 ) 



(64) 



with D = [(1 + n 2 )A]- 1 and AA = A+ - A". De- 
pending on the values of the parameters n\, n 2 , and 
n 3 the total energy density (fBTij) behaves asymptotically 
as ph — > a 3A /( 1+ ™ 2 ) in the limit of large scale fac- 
tors, meaning that a t-2(i+« 2 )/3A± _ For c = (+) 
or 6 =0 (-), the model includes the exact power-law 
expansions = £- 2 (i+" 2 )/3A _ j n particular, when 
(l + n 2 )/3A ± > 0, we have a final phantom phase. In the 
special cases that A ± = and A T < 0, namely n\ > 
and 7J3 = or ri\ < and 713 = (see Eq. (|64I0 . we 
have a final de Sitter stage with an effective cosmologi- 
cal constant given by the limit ph — > h-eff = 6 1 /( 1+ ™ 2 ' 
or ph — > A e // = c 1 ^ 1 " 1- " 2 ). These models include the 
modified Chaplygin gas introduced in Ref. where it 
was proposed the equation of state p = Ap — B / p n , with 
n > 1, and the parameters A and -B were constrained 
to be positive. When both, A + 7^ and A~ 7^ 0, the 
equation of state (|63[) contains those which characterize 
various of the variable modified Chaplygin gas models 
investigated in Refs. [l]| -0. 



B. Nonlinear examples 

Here we examine the case where Qh takes the form 

PcPx 



Qho = O-lc 



(65) 



with a constant, see [6j and references therein. From 
Eqs. (|58|) and (|65p we identify the three coefficients n\ = 
(7c + 7x)(l + n 2 ), n 2 = aj c /A and n 3 = ~f c j x (l + n 2 ). 
Then, the characteristic roots read 



-7x(l + n 2 ), A- = -7c(l + n 2 ). (66) 



Introducing these roots into the general solution and 
the effective equation of state they reduce to 



ca -3 7c (l+ri 2 ) +&a -3 7 x(l+«2) 



l/(l+n 2 ) 



p= (7 C - l)p- 6 A 



a -3 7a ,(l+n 2 ) 



(67) 



(68) 



This effective equation of state can be identified with 
those which were used to build variable modified Chap- 
lygin gas models [l8|-(23l|. 



1. Modified Chaplygin gas 

We make an adequate selection of the parameters so 
that the interaction ([6"5")) be focused on a dark energy 
component described by some kind of vacuum energy 
density i.e., ~f x = 0. Then, n 3 = j c Jx(l + n 2 )=0, A = j c , 
n 2 = a =/= — 1, A + = 0, and the second term in Eq. 
([67| becomes constant. Assuming that the dark mat- 
ter component is nearly pressureless, we may associate 
it to a barotropic fluid with a free constant parameter 
7c = 70 ~ 1- Then the energy density (|67]) abbreviates 
to 



P = 



B /ao\ 3 T ( 1+Q! ) 
7o v a ' 



1/1+a 



(69) 



where B is a constant. Hence, by replacing this energy 
density in Eq. (|9|) , we obtain the equation of state of the 
one effective fluid 



P = (70 - l)p~ : 7 ^- 
P 



(70) 



It characterizes several unified cosmologies implemented 
with Chaplygin gases as the generalized, extended, mod- 
ified and enlarged ones 0-[I3- 

Now, we express the Chaplygin gases with equation 
of state ([70)1 as an interacting two-fluid model where the 
energy exchange is produced by the interaction term (|65[) . 
To this end, we insert the equation of state ([70j) into Eq. 
([5]) to find the dark matter and dark energy densities 



PcCh 



B 

7o/0° 



PxCh 



B 

loP° 



rc h = -l + ^P a+1 . 



(71) 



(72) 



On the other hand, by combining the Eq. (fTTj) 
and (|72[) . we can express the quantities characterizing the 
interacting two-fluid model in terms of the scale factor. 

For expanding universes and 70(1 + a) > 0, the total 
energy density (f6"9"| p a+1 —> B/jo and the ratio has a 
vanishing limit rch 0. 



2. The 70 = case 

When 70 = 0, the effective one-fluid model described 
by the expressions (|69)) and ([70)) is not valid. However 
this case can studied from the nonlinear interaction term 



Q, 



a(k 4 - l)p V 2 + "idxP + he + lx)p' + k 4 p" 



(73) 

Now, the source equation (|59j) reduces to x" = and one 
finds that the total energy density p = xVU+c) 



Po 



±1 + b In 



1/(1+") 



(74) 



8 



has a logarithmic dependence with the scale factor, while 
the equation of state of the effective one-fluid model is 



P 



-P 



1 + Q 



(75) 



with po a constant. The last equation of state also can 
be seen as a generalization of the polytropic equation of 
state p = Kp" <p , where K is a constant and 7 P is the 
polytropic index. The scale factor is determined by inte- 
grating the Friedmann equation for the source (|74[) 



1 

a = a exp — 
36 



±1 



2(1 + a) 



(76) 

where we have set to = 0. Note that the quantities ([71)1 - 
(|76p are independent of k±, then the interaction term p" 
dos not contribute to the evolution of this unified model. 



C. The "relaxed Chaplygin gas model" 

Here we are going to consider the "inhomogeneous non- 
linear interaction" (Qi), 

fcip-y 2 + k 2P + k 3P ' + k AP " + /(77)p fei /(i-M 



A 



Q l = 

(77) 

which has been obtained from the Eq. ([57]) by choos- 
ing v = —ri2 = ki/(l — k±). The source equation ([59]) 
becomes inhomogeneous and its general solution is given 
by 



ca 3X + ba 3X+ +x r 



l/(l+n 2 ) 



while the equation of state ((9]) takes the form 



Pi = 



1 + n 2 



bAX 



,3a" 1 



A- 



Pi 



II + n 2 )p r i 



(78) 



(79) 



(i) For f = fa = const the particular solution x p = 
n 4/o/^3 is constant and the total energy density (|78[) 
along with the effective equation of state ([79)) describe a 
"double unified model" , in a sense that initially the uni- 
verse is dominated by the "two" terms inside the bracket 
of pi w (c« 3r + ba 3X+ ) 1 ^ 1+n2 \ These terms can be 
seen as a "nonlinear mixture of two fluids". But at 
late times the universe is dominated by a vacuum en- 
ergy pi w Xp 1+n2 ^ and has a de Sitter expansion. So, 
this (i) case includes various generalizations of the mod- 
ified Chaplygin gas model investigated in Refs. fl8j|-[23l| . 
When one of the constants c or b vanishes the effective 
equation of state ([79"]) turns in Eq. ([70)) and the double 
unified model produces different versions of the Chaply- 
gin gas [11- [l3| . 

(ii) For f(rj) 7^ const and 77 = In a, the particular solu- 
tion x p and consequently the numerator of the last term 



in the RHS of the effective equation of state ([79]) be- 
come arbitrary functions of the scale factor. Then, a 
fluid obeying the Eq. ([79")) defines a "relaxed Chaplygin 
gas model" . 



V. CONCLUSIONS 

We have investigated interacting dark sector models 
with energy transfer and introduced an effective one-fluid 
description with an effective equation of state. The two 
coupled equations describing the interacting model have 
been combined to obtain the fundamental second order 
differential equation for the total energy density, "the 
source equation" . We have assumed a separable interac- 
tion Q — p<9(7,7')j which includes a large set of cases 
investigated in the literature, and examined the stability 
condition for the power-law solutions. 

We have presented the "linear interaction" Ql = 
cip c + C2p x + C3p' c + C4p' x + c 5 p+c 6 p' + c 7 p" , and reduced it 
to a linear combination of the basis elements p, p' and p" , 
so that Ql = b\p + b 2 p' + b 3 p" . As far as we know, this 
interaction has not been investigated in the literature. 
We have found the stationary solution for the effective 
barotropic index, 7 S and "f + , and imposed the stability 
condition, being the power-law solution a = t 2 / 3ls an at- 
tractor. Interestingly, the existence of the j s solution is 
linked to the fulfillment of the requirement Ql{1s) < on 
the attractor, indicating that the energy is being trans- 
fered from the dark energy to the dark matter. We have 
considered several particular examples and observed that 
Ql <x Px-, Ql oc p and Ql oc p' are satisfactory couplings 
because in each case, the ratio r s = p cs / 1 p xs is an attrac- 
tor. These simple models may alleviate the problem of 
the coincidence. Although, a coupling proportional to 
the dark matter energy density does not lead to stable 
solutions however, it can work when it is combined with 
the remaining terms of Ql- Finally, we have obtained 
the exact scale factor and the effective barotropic index 
in implicit form. 

We have generalized the above coupling by adding a 
constant to the Ql and introduced the "general linear 
interaction", Q gL = A eff (-f c j x - &]A)/A + Q L - These 
models lead to several alternatives to the ACDM model. 

We have presented a class on "nonlinear interaction" 
QnL = [fcip-V 2 + fep + fc3P' + fc4/o" + /(r7)/ l/(1 - fe4) ]/A. 
Although the source equation becomes a nonlinear dif- 
ferential equation, we have linearized and solved it. In 
general, we have found that the equation of state of the 
effective one-fluid model depends explicitly on the scale 
factor and the universe evolves to a power-law scenario 
for large cosmological times. However for interactions 
having the form Q ti l = ori c p c p x j p with a a constant, we 
have shown that the effective equation of state becomes 
that of the Chaplygin gas when the dark energy compo- 
nent is described by some kind of vacuum energy density 
i.e., 72; = 0. These unified model has been expressed as 
an interacting one. Also, we have examined the partic- 
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ular nonlinear interaction leading to the source equation 
x" =0. Here we have shown that the effective equation 
of state can be seen as a generalization of the polytropic 
equation of state. 

Generically, when there are no restrictions on Q n L, the 
equation of state of the effective one-fluid model defines 
what we have called the "relaxed Chaplygin gas model" . 
It contains various generalizations of the Chaplygin gas, 
including the variable modified Chaplygin gas model with 
equation of state p = Ap + B{a)/ p a . 
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